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Abstract 

In this paper we describe a new release of the PRESHOWER program, a 
tool for Monte Carlo simulation of propagation of ultra-high energy photons 
in the magnetic field of the Earth. The PRESHOWER program is designed 
to calculate magnetic pair production and bremsstrahlung and should be 
used together with other programs to simulate extensive air showers induced 
by photons. The main new features of the PRESHOWER code include a 
much faster algorithm applied in the procedures of simulating the processes 
of gamma conversion and bremsstrahlung, update of the geomagnetic field 
model, and a minor correction. The new simulation procedure increases the 
flexibility of the code so that it can also be applied to other magnetic field 
configurations such as, for example, encountered in the vicinity of the sun or 
neutron stars. 
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PRESHOWER 2.0 

ADWG_v2_0 

[http: / / cp c . cs . qub .ac.uk/ summaries / 
ADWG_v2_0.html 

CPC Program Library, Queen's Uni- 
versity, Belfast, N. Ireland 
Standard CPC licence. 



1. Program Summary 

Program title: 
Catalog identifier: 
Program summary URL: 

Program obtainable from: 

Licensing provisions: 

Programming language: 

Computer(s) for which the program has 

been designed: 

Operating system(s) for which the pro- 
gram has been designed: 
RAM required to execute with typical 
data: 

CPC Library Classification: 
External routines/libraries used: 
Catalog identifier of previous version: 
Journal Reference of previous version: 

Does the new version supercede the 
previous version?: 
Nature of problem: 

Solution method: 



Reasons for the new version: 
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[http: / / cp c . cs . qub .ac.uk/licence/ licence . html 
C, FORTRAN 77 
Intel-Pentium based PC 

Linux or Unix 

< 100 kB 

1.1. 

IGRF [11,3, DBSKA Q, ran2 H 
ADWG_vl_0 

Computer Physics Communications 

173 (2005) 71-90 

yes 

Simulation of a cascade of particles 
initiated by UHE photon in magnetic 
field. 

The primary photon is tracked until its 
conversion into an e^e~ pair. If conver- 
sion occurs each individual particle in 
the resultant preshower is checked for 
either bremsstrahlung radiation (elec- 
trons) or secondary gamma conversion 
(photons) . 

1) Slow and outdated algorithm in the 
old version (a significant speed up is 
possible); 2) Extension of the program 
to allow simulations also for extrater- 
restrial magnetic field configurations 
(e.g. neutron stars) and very long path 
lengths. 



Summary of revisions: 



Restrictions: 



Running time: 



A veto algorithm was introduced in the 
gamma conversion and brcmsstrahlung 
tracking procedures. The length of the 
tracking step is now variable along the 
track and depends on the probability 
of the process expected to occur. The 
new algorithm reduces significantly the 
number of tracking steps and speeds 
up the execution of the program. The 
geomagnetic field model has been up- 
dated to IGRF-11, allowing for interpo- 
lations up to the year 2015. Numerical 
Recipes procedures to calculate modi- 
fied Bessel functions have been replaced 
with an open source CERN routine DB- 
SKA. One minor bug has been fixed. 
Gamma conversion into particles other 
than an electron pair is not considered. 
Spatial structure of the cascade is ne- 
glected. 

100 preshower events with primary en- 
ergy 10^° eV require a 2.66 GHz CPU 
time of about 200 sec; at the energy of 
10^1 eV, 600 sec. 
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2. Introduction 



Identifying and understanding the sources of cosmic rays with energies up 
to 10^° eV is one of the most important questions in astroparticle physics (see, 
for example, js], 0, 0]). Knowing the fraction of photons in the flux of ultra- 
high energy cosmic rays is of particular importance as photons are unique 
messengers of particular source processes (acceleration vs. decay of super- 
heavy particles or other objects). They are also produced in interactions 
of charged cosmic ray nuclei of the highest energies with cosmic microwave 
background radiation, known as the Greisen-Zatsepin-Kuzmin (GZK) effect. 
If the energy of the cosmic ray particles exceeds the GZK energy threshold at 
the source, ultra-high energy photons are produced due to well-understood 
hadronic interactions with microwave photons and can be detected at Earth 
as a unique propagation signature. So far only upper limits to the photon 
flux at ultra-high energy exist, which have led to severe constraints on models 
for UHECR sources. However, depending on the primary cosmic ray com- 
position, the sensitivity of the latest generation of cosmic ray detectors, i.e. 
the Pierre Auger Observatory [sj and the Telescope Array jl], should allow 
detection of GZK photons for the first time. 

The simulation of the propagation of photons before they reach the Earth's 
atmosphere is important because of the preshower effect 10| that may occur 
when a photon traverses a region where the geomagnetic field component 
transverse to the photon trajectory is particularly strong. As described e.g. 



in Refs. [111 . |12| , high energy photons in the presence of a magnetic field may 
convert into e^e~ pairs and the newly created leptons emit bremsstrahlung 
photons, which again may convert into e"'"e~ if their energies are high enough. 
As a result of these interactions, instead of a single high energy photon, a 
shower of particles of lower energies, the so-called preshower, reaches the 
atmosphere. The occurrence of the preshower effect has a large impact on 
the subsequent extensive air shower development and changes the predicted 
shower observables. 

In 2005 the program PRESHOWER for simulating the showering of 
photons in the Earth's magnetic field was released. It was shown that this 
initial version of the Monte Carlo code is in good general agreement with 



previous studies [10|, ll3|, [ij, llSl, lla, ll7|, ll8| . In this paper we describe the 



changes of the PRESHOWER code relative to the initial version 1.0 [11 
The main feature of the new release is a much faster algorithm for calcu- 
lating the distance at which a preshower interaction (gamma conversion or 
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bremsstrahlung) occurs. In version 1.0, the calculations were done in con- 
stant steps along the particle trajectory and the step size was optimized for 
all possible trajectories. The constant step practically disabled studying the 
preshower effect along paths longer than several tens of thousands kilometers. 
Now, in version 2.0, the distance to the next interaction point is computed 
with an efficient veto algorithm, decreasing significantly the number of com- 
puting steps. The new method, being independent of the trajectory length, 
allows computations for arbitrarily long photon paths, e.g. simulations of 
preshower creation in the vicinity of a neutron star or active galactic nu- 
cleus. The new algorithm is described in detail in Section |3] of this article. 
Other important changes, i.e. the update of the geomagnetic field model and 
a code correction in version 1.0 are discussed in Section SJ The results of 
testing PRESHOWER 2.0 are presented in Section [5] and conclusions are 
given in Section |6l 

Following Ref . ll| , all the results presented in the following are obtained 
for the magnetic conditions of the Pierre Auger Observatory in Malargiie, 
Argentina (35.2°S, 69.2°W). The shower trajectories are given in the local 
frame where the azimuth increases in the counter-clockwise direction and 
= 0° refers to a shower coming from the geographical North. 



3. The new sampling algorithm 

In PRESHOWER 1.0 the effect of precascading is simulated following the 
particle trajectories with a fixed step size. In each step the probability of 
conversion into e^e~ is calculated for photons and the probability of emitting 
a bremsstrahlung photon is computed for electrons. The step size has to 
be optimized for all possible trajectories and magnetic field configurations 
encountered along the particle trajectories. It has been found out that a 
step size of 10 km works well until the primary photon conversion and then 
the emission of bremsstrahlung photons as well as conversions of secondary 
photons are simulated in steps of 1 km. In this algorithm a typical simulation 
run consists of several thousands of steps. 

The number of simulation steps and hence the computing time can be 
significantly reduced by using a veto algorithm. The algorithm used in the 
following is commonly applied in physical situations where a probability of 
occurring of a certain process varies within a given spatial or temporal inter- 
val (see, e.g. Ref. [loj). The location or time of the occurrence of a physical 
process studied is found in few approximating jumps. An example of a pro- 
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cess that can be treated with the algorithm is a radioactive decay considered 
in a certain interval of time. Photon conversion and bremsstrahlung proba- 
bilities to be found along spatial trajectories can also be computed with this 
veto algorithm. 

3.1. General description 

The theory behind the veto algorithm is based on a process of discrete 
events described by 

^| = -/W^W. (1) 

Then the probability cIPa of occurrence of process A in the time window 
t . . .t -\- dt is given by 

Together with the probability of not having an occurrence of process A in 
the time from to to t 

P (3) 
^"""-^ ~ N{to) ^ ^ 

one obtains for the probability dP for having an occurrence of A in the time 
window t . . .t + di, provided that this process did not occur at an earlier time 
t' with to <t' <t 

dP = Pno-A dPA = fit) dt = fit) exp |- ^ fit') rft j dt . (4) 

If an analytic solution can be found for the integral of fit) 

Fit)= f fit')dt' (5) 

one can sample the time t of the next occurrence of A after the previous 
occurrence at time to using the inversion method 



f dP = exp{F(i)} = e, t^ F-\\ni), 



(6) 



with ^ being a random number uniformly distributed in (0, 1]. 

If Fit) cannot be found or the inverse of it computed sufficiently easily 
one can use a function git) such that ^t > : git) > fit) and apply the 
following procedure 
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1. set the initial conditions: i = 0, to = 0; 

2. i = i + l; 

3. get a random number G (0, 1); 

4. compute U = G~^{G{ti^i) - In^j), U > 

5. get another random number C,i G (0, 1); 

6. if < f{ti)/g{ti) then ti is the wanted result, i.e. the moment when 
process A occurred, otherwise one has to go back to step 2 or, if the 
end of the interval is reached, end the procedure without occurrence of 
process A. 

The algorithm described above was mathematically proven to reproduce 



exactly the expected distributions [20 



3.2. Implementation in PRESHOWER 2.0 

Following the general scheme described above, the application of the veto 
algorithm to simulations of the preshower effect is straightforward. Instead of 
the time variable t, the distance r along the preshower trajectory is used. We 
consider two processes in parallel over an interval starting at r start and ending 
at Tend-, namely gamma conversion and bremsstrahlung of electrons. Follow- 
ing the physics notation introduced in Ref. [ll| (see also the Appendices 
A and B for all the required physics formulas and symbols) the probability 
functions are defined as 

Pconvir) = a{xir)) (7) 
(see Eqs. IA.ltlA.5|) for gamma conversion and 



P,rem{r)^ ! I {B ^{v) , E M)'^ (8) 







hv 



(see Eqs. IB.ltlB.7p for magnetic bremsstrahlung. The function f{t) is then 
replaced by Pconv{r) or pbremij), depending on the process to be simulated. 
Since finding the antiderivatives of Pconvij') or Pbrem{r) is not straightforward, 
simple functions limiting Pconvi^) and Pbremij) are used. We define 

gconv{r) = p'^^Z = const, Vr G {Vstart, Tend) ■ PTcZ > Pconv{r), 
9brem{r) = pZem = COHSt, Vr G (r star t , r end) ■ PZem > Pbremir), 



which replace g{t) and for which the antiderivatives are 

G conv 

^hrem\i } Fbrem ' • 



(10) 
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With these substitutions the algorithm of Sec. 13. H is apphed in PRESHOWER 
2.0 for the two interaction processes. 

The determination of V^^vibrem crucial for computing time optimiza- 
tion. Too large p]^onv/brem incrcases the total number of steps to be executed 
in the procedure. 

3.3. Determination of p^^^^^^^^^ 

The functions Pconv/bremij) depend on -B±(r), which is computed with a 
numerical model. Hence finding the absolute maxima p^"^ and p^em done 
numerically. Moreover, through the dependence on -B±(r), both Pconvij") 
and Pbremi'^) depend on the primary arrival direction and the geographical 
location of the observatory. As can be seen in Figs. [T] and [21 the values of 
p^^y and p^r'em be significantly different for various arrival directions. 

This indicates that the computation of p^^^ and p|^em should be per- 
formed for each trajectory separately, otherwise one would have to apply up- 
per limits of these values which would be universal but excessively large for 
most directions. Accepting the excessive values of V^^v/hrem would increase 
enormously the number of steps in the veto algorithm and might result in an 
unacceptable increase of the computing time. 

Typically, the functions Pconvir) and Pbremij) reach their global maximum 
at the top of the atmosphere, i.e. at the end of preshower simulations, as- 
sumed here to be at the altitude of 112 km. This is the point closest to the 
Earth's surface and Bx{r) typically reaches the maximum value. However for 
certain classes of trajectories -B±(r) might start to decrease with approach- 
ing the geomagnetic field source when the trajectory direction approaches 
a tangent to the local field lines. If this decrease happens to be close to 
the Earth surface, preshower particles are exposed to the maximum -B_L(r) 
somewhere before reaching the atmosphere. Examples of p^^^ and pJ^eJ^ with 
local extrema are plotted in Figs. [1] and [2l The positions of the minima of 
Vconv/hremij') are closcly correlated with the minima of B^{r). 

In case of Pconv{f)i it has been checked that its global maximum is well 
reproduced by computing the function value along the trajectory in a simple 
loop with steps of 1000 km. This procedure is performed for the primary 
photon energy and the trajectory of interest. It has been checked that Pconvir) 
decreases with energy, so p^^^ found for the primary photon energy will work 
also for secondary photons of lower energies. In case of electrons, Pbrem{r) 
may increase with decreasing electron energy, so one has to compute p^em ^^r 
energies within the entire energy range of the simulated particles. Here p^g^ 
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Figure 1: Examples of Pconv functions along different trajectories at the location of the 
Pierre Auger Observatory in Malargiie (Argentina) . A minimum value for one of the curves 
is related to the small value of B± for this specific arrival direction and altitude. See text 
for further details. 
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altitude a.s.l. [km] 

Figure 2: Examples of pbrem functions along different trajectories at the location of 
the Pierre Auger Observatory in Malargiic (Argentina). A minimum value for one of the 
curves is related to the small value of B± for this specific arrival direction and altitude. 
See text for further details. 
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is computed in two loops. The external loop along the trajectory is done in 
steps of 1000 km and the internal loop over energies decreases the energy by 
one decade in each step. The steps in both procedures can be adjusted by 
the user if necessary. There is also an alert in the program that gets triggered 
when the actual values of Pconv/bremir) happen to exceed p'^^l/hrem- 

The above method of finding the absolute maxima of Pconv/brem{r) is fast 
and efficient. However, it is optimized only for specific simulation conditions: 
preshowering in the geomagnetic field. In other environments, involving more 
irregular shapes of B{r), one has to reconsider the procedure of finding the 
absolute maxima of Pconv/bremij')- 

4. Other modifications and corrections 

Other modifications and changes applied in the new release of the PRESHOWER 
program are briefly listed below. 

1. The IGRF geomagnetic field model has been updated to the year 2010 
and the most recent IGRF- 11 coefficients have been applied (Ref. (2l|). 
In the updated model, the highest order of spherical harmonics has been 
increased from 10 to 13. The geomagnetic field can be extrapolated up 
to the year 2015 with the new model. The differences between the field 
strength and direction in these two models are not larger than 0.001%. 

2. Procedures to calculate modified Bessel functions have been replaced 
with an open source CERN routine DBSKA. 

3. A minor problem has been found in the auxiliary function kappa(x) 
used for calculation of bremsstrahlung probability. The interpolation 
performed in this function failed for the rare case of x = 10.0. This 
happened because of a faulty definition of the last interval where the 
interpolation was done. As a result of this bug the input value x = 10.0 
was excluded from the computations. This bug has been fixed in the 
new release of the program. 

4. Since for some cases the number of preshower particles can be very 
large, the size of the array part_out [50000] [8] , which stores the out- 
put particle data, has been increased from 50000 to 100000 entries. 

5. The code of the program was reorganized and more clearly structured. 
The main change here was moving the auxiliary functions and routines 
to a separate file. 

A list of the new and modified files with basic explanations can be found in 
the Appendix C. 
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5. Validation of the new version 

The new release has been intensively tested. Below we show some exam- 
ples to illustrate the performance. 

In Fig.[3]a comparison of conversion probability obtained with PRESHOWER 
1.0 and PRESHOWER 2.0 is presented for different arrival directions and 
primary energy of 7x 10^^ eV. The lines represent conversion probabilities ob- 
tained by numerical integrations of the expression (IA.5|) along trajectories a 
nd with steps as given in PRESHOWER 1.0. The points are plotted to show 
fractions of events with gamma conversion simulated by PRESHOWER 2.0. 
Each fraction was calculated after 10,000 simulation runs. Simulations of 
gamma conversion probabilities for other primary energies has also been 
checked and in all cases an excellent agreement between the results of the 
two PRESHOWER versions has been found. 

A cross-check of the procedures responsible for simulation of bremsstrahlung 
is shown in Fig. HI in which the energy distribution of secondary particles for 
a primary photon of 10^° eV and an arrival direction along which the trans- 
verse component of the geomagnetic field is particularly strong ( "strong field 
direction" ) are compared. Plotted are the summed distributions of energies 
of secondary photons and electrons together with the relevant histograms 
weighted by the energies. The summations are done for 500 simulation runs. 
The results obtained with the two program versions are in very good agree- 
ment. 

Further tests for the same set of simulations are shown in Figs. |5] and [61 
These are the number of preshower particles and the total energy carried by 
the preshower electrons. Both observables are calculated at the top of the 
atmosphere and both are plotted versus the altitude of primary photon con- 
version. In both figures a comparison is made between the results obtained 
with PRESHOWER 1.0 and PRESHOWER 2.0. Again, the agreement be- 
tween the results of the two PRESHOWER versions is very good. 

One of the main aims of the new release of PRESHOWER was to reduce 
the computing time. The results of the CPU time comparison are summa- 
rized in Table [H The computing time is more than a factor 5 shorter in case 
of simulations with PRESHOWER 2.0 than in the case of PRESHOWER 1.0. 
This reduction is seen both in computation of gamma conversion (speed up 
by nearly factor 10) and in more time consuming bremsstrahlung routines. 
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50 100 150 200 250 300 350 



azimuth [deg] 

Figure 3: Total probability of 7 conversion for the primary energy of 7 x 10^^ eV for dif- 
ferent arrival directions as computed by PRESHOWER 1.0 (lines) and PRESHOWER 2.0 
(points). PRESHOWER 1.0 values were obtained by numerically integrating the conver- 
sion probability in the loop over trajectory. The points represent the fractions of events 
with gamma conversion simulated by PRESHOWER 2.0 with the new veto algorithm. 
Each fraction is the; average for 10000 primary photons. Computations have been done 
for magnetic conditions at the Pierre Auger Observatory in Argentina. The azimuth 0° 
refers to showers arriving from the geographic North. 
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Figure 4: Energy distribution of photons (top left) and electrons (bottom left) in 500 
preshowers initiated by 10^° eV photons arriving at the Pierre Auger Observatory in 
Argentina from the strong field direction. The spectra weighted by energy are plotted to 
the right. The dashed histograms were obtained with PRESHOWER 1.0 and the solid 
ones with PRESHOWER 2.0. 
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4500 



first conversion altitude a.s.l. [km] 

Figure 5: Number of particles in the preshower for different altitudes of the first 7 conver- 
sion simulated with PRESHOWER 1.0 and PRESHOWER 2.0. Plotted are the preshowers 
initiated by 10^" eV photons arriving from the strong field direction. The two points some- 
what higher than those of the general trend are cases where one of the bremsstrahlung 
photons again converted in the magnetic field to produce an electron-positron pair which 
emitted the additional photons. See also Fig. [S] 

Table 1: The preshower simulation times in PRESHOWER 1.0 (old) and PRESHOWER 
2.0 (new) for selected arrival directions and primary energies. The arrival directions are 
selected to represent a typical variation of B±: 6 = 70° and = 0° for a "weak" B±, 
0^0° for a "medium" B_l and 9 = 60° and (p = 177° for a "strong" Bj^. 



Eo [eV] 


direction 


fraction of 


time old 


time new 






converted 


[sec] 


[sec] 


7x10^9 


6 = 0° 


0/1000 


79 


7 


7x10^9 


e = 70°, = 0° 


0/1000 


76 


8 


1020 


e = 60°, = 177° 


92/100 


1195 


209 
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Figure 6: Energy carried by preshower electrons at the top of the atmosphere vs. the 
ahitude of the first 7 conversion for a primary photon energy of lO^'^ eV in the strong field 
direction. The two points in excess of the general trend are two rare cases where the first 
bremsstrahlung photon converted again into an electron-positron pair which increased the 
total energy carried by leptons. See also Fig. [5j 
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6. Summary 



The program PRESHOWER is a tool designed for simulating magneti- 
cally induced particle cascades due to ultra-high energy photons. It can be 
linked with air shower simulation packages such as CORSIKA 22(| to calcu- 
late complete photon-induced particle cascades as they are searched for with 
cosmic ray observatories. 

A new version of the PRESHOWER program, version 2.0, has been re- 
leased and its features are presented in this article. An efficient veto algo- 
rithm has been introduced to sample the locations of individual interaction 
processes. Other modifications include the update of the geomagnetic field 
model, correcting a rare exception, and reorganizing the program code. The 
results obtained with the new release of PRESHOWER agree very well with 
those calculated with the previous version. 

The new algorithm not only speeds up the program by more than a factor 
5, but also allows additional applications due to the increased fiexibility of 
the sampling of interaction points. For example, the preshower effect can 
now be studied not only in the geomagnetic field, but also close to extended 
astrophysical objects like neutron stars and active galactic nuclei. An appli- 
cation of PRESHOWER 2.0 in the conditions other than the geomagnetic 
field require only small changes in the program. The magnetic field model has 
to be replaced and start and end points of simulations have to be adequately 
adjusted. 
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Appendix A. Magnetic pair production: 7 — )■ e+e 

The number of pairs created by a high-energy photon in the presence of a 
magnetic field per path length dr can be expressed in terms of the attenuation 
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coefficient a{x) |23|: 

ripairs = nphotonsi'^ " exp[-a{x)dr]} , (A.l) 

where 

a{x) = 0.5{aemmec/h)iBjB,r)T{x) (A.2) 

with aem being the fine structure constant, x = 0.5{hh'/meC^){Bj_/Bcr), B± 
is the magnetic field component transverse to the direction of the photon's 
motion, Bcr = m^c^/eh = 4.414 x 10^^ G and T{x) is the magnetic pair 
production function. T(x) can be well approximated by: 

T(x) = 0.16x-'K\/,{^), (A.3) 

where -ft'i/s is the modified Bessel function of order 1/3. For small or large 
arguments T{x) can be approximated by 



nx) ^ ^ (a.4) 



0.46exp(-A), 

0.60x-^/^ X > 1- 

We use Eq. (lA.ip to calculate the probability of 7 conversion over a small 
path length dr: 



Pconv{r) = 1 - exp[-a{x{r))dr] ^ a{x{r))dr. (A.5) 



Appendix B. Magnetic bremsstrahlung 

After photon conversion, the electron-positron pair is propagated. The 



energy distribution in an e~^e pair is computed according to Ref. [24 



da{e,x) _ aemmecB± 3^/^ [2 + e{l - e)] 
de hBcr 97rx e(l — e) 

where e denotes the fractional energy of an electron and the other symbols 
were explained in the previous chapter. The probability of asymmetric energy 
partition grows with the primary photon energy and with the magnetic field. 
Beginning from x > 10, the asymmetric energy partition is even more favored 
than the symmetric one. 

Electrons traveling at relativistic speeds in the presence of a magnetic 
field emit bremsstrahlung radiation (synchrotron radiation). For electron 



?,Xe{l-e) 



(B.l) 
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energies E ^ mgC^ and for B± <^ B^r, the spectral distribution of radiated 
energy is given in Ref. [251] : 

where ^ = {3/2){B±/Bcr){E/meC^), E and me are electron initial energy and 
rest mass respectively, K^/^ and K2/3 are modified Bessel functions, and y is 
related to the emitted photon energy hu by 

The total energy emitted per unit distance is (in CGS units) 

W = \rlBl(^X r f{y)dy (B.4) 

with To being the classical electron radius. For our purposes we use the 
spectral distribution of radiated energy defined as 

, , hvdN ,^ , 

where dN is the number of photons with energy between hv and hv + d{hv) 
emitted over a distance dx. From Eqs. (]B.2p . (IB.Sp . (lB.4p . and (IB.SP we geto 

HBM - \rlBl . (B.6) 

Provided dx is small enough, dN can be interpreted as a probability of emit- 
ting a photon of energy between hv and hv^dihv) by an electron of energy E 
over a distance dx. In our simulations we use a small step size oi dx =\ km. 
The total probability of emitting a photon in step dx can then be written as 

Pbrem{B±,E,hu,dx) = [ dN = dx [ I{B±,E,hu)^^^l^ . (B.7) 

J Jo hi' 

The energy of the emitted photon is simulated according to the probability 
density distribution dN/d{hv) obtained from Eq. IB. 71 



Expression (jB.6|) . valid for all values of hv, is equivalent to Eq. (2.5a) in Ref. 23| . 



A simplified form of distribution 
be used only for hv <C E. 



Re]) is given by Eq. (2.10) in Ref. 



however it can 
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Appendix C. Details on the new and modified files included in the 
PRESHOWER 2.0 package 



The files included in PRESHOWER 2.0 package but not existing in the 
previous release are : 

IGRF-11 . f An external routine generating the geomagnetic field components 
according to the lGRF-11 model [2|. This file replaces igrf .f of the 
previous release of PRESHOWER. 

cernbess . f External procedures calculating the sequence of modified Bessel 
functions [3]. These open source procedures replace previously used 
functions from Numerical Recipes. 

veto . c This file contains functions and procedures called by the veto algo- 
rithm. 

veto .h The header file for veto . c. 

utils.c This file contains auxiliary functions and procedures used within 
the program. In the previous version of PRESHOWER the auxiliary 
functions were placed in preshw.c, now it is more convenient to have 
them in a separate file. 

utils . h The header file for utils . c. 

The list of files which existed in the previous release of PRESHOWER 
but have been modified for PRESHOWER 2.0 include: 

preshw.c contains the main procedure preshw_veto generating preshowers 
with the veto algorithm, 

prog.c reads input parameters and calls preshw_veto, 

Makefile modified to account for the new files. 
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